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ABSTRACT 


This report exanlnes some aspects of the lack of close correlation between 
the predictions of analytical nodeling of dynanic structures and the results of 
vibration tests on such structures, and suggests ways in which the correlation 
may be improved. 
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SECTION I 
INTHOnrCTION 


Desiftning a structure, to survive a prescribed dynamic environment is 
most often performed nowadays by usin;; an analytical model ot t!ie structure- 
Since many structures will not be subiected to their design environment prior 
to commissioning, it is very important, therefore, that the analytical code 
mimic the behavior of the physical system with reasonable accuracy. KiUle 
modern analytical techniques have the capability ot modeling a physical system 
to any desired degree of accuracy, financial or time conslderat ions may pre- 
clude the use of liighly accurate analytical models of the system. Clearly, if 
the cost of analytically modeling the system becomes comparable with tl:e cost 
of building and testing tlie physical system, then it may well be that cut-and- 
try methods are more cost effective than analytical methods. While the air- 
craft and aerospace industries have used analytical modeling for decades, 
until recently the automotive Industry found that, due to the complex double 
curvature of the automobile body, it vras more cost efficient to use cut-and— trv' 
methods rather than analytical modeling- To compromise between cut-and-try 
methods and the use of highly accurate, but very expensive, analytical models, 
engineers arc frequently willing to accept a fairly crude analytical model of 
the desired structure and to resort to a limited program of testing to 
"qualify” the analytical model. Since the engineer, unllhe the mathemat Ician 
or scientist, must always balance rigor against cost, an important question 
is, "Given that the d\'nnnic environment is known with only limited accuracy, 
how accurate must the analytical model bo to obtain ’adequate* predictions of 
the dynamic response of the structure?** This report examines a number of 
aspects of the problem of trying to correlate the results of dynamic testing of 
a structure with the analytical predictions based on rather crude modeling. 



SECTIOK II 

A^*ALYTIC/VL MODELING AND PREDICTION OF RESPONSE 

Almost all structural systems are distributed parameter (continuous) 
systems; this is particularly true in the case of aerospace vehicles where 
the desire to nlninize veight results in a design with siass and stiffness 
distributed throughout the system. The complex geometry and boundary condi- 
tions in space vehicles seldom permit exact solutions of the partial dif- 
ferential equations describing the dynamical behaviors of the vehicle. For 
this reason, finite differences, finite elersent, Raylcigh-Ritz, or Calerkin 
techniques are normally used to discretize the system and reduce the problem 
to that of a lumped parameter system. VHille these techniques differ in 
detail, they all have the same general properties. They attempt to approxi- 
mate a space-continuous system by a discrete system having a finite number of 
degrees of freedom, N. A common feature of such discrete approximating 
schemes is that only the first M modes, M = 1/3 N, have a reasonable chance of 
accurately modeling the first M modes of the continuous system. Katurally, 
the larger N is, the larger H nay be, and the better the degree of approxima- 
tion in the lover modes. The higher discrete modes are, in general, poor 
approximations of the continuous modes, even when N is large. 

Since the stress in a continuous structure depends on the spatial deriv- 
atives of the deformations, the order N of the discrete approximating system 
should be large to obtain accurate approximations of the spatial derivatives; 
however, in many important problems in structural dynamics, the number of 
dynamically active nodes, «, is much smaller than N. For this reason the 
analyst will frequently reduce the size of his model to core closely corre- 
spond to the number of active modes. This practice, unfortunately, reduces 
the accuracy with which stresses and forces in the structure cay be deter- 
mined, particularly in the case of transient notions. 


A. STRUCTURAL DYNAMICS 
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daaping and admits classical nonaal nodes. It will be shown that, if the 
danping is snail and the eigenvalues well separated, this is a reasonably good 
approxinjation, at least in the lover inodes. 

1. Continuous Systems 

To illustrate the techniques, let us restrict our attention to a rela- 
tively simple continuous structure such as a beam, a plate, or a shell that 
can be described by Equation 2-1: 


+ L^u = f(x, t) on D 
u(x, 0) = 0) = 0 


\ 


> 


(2-1) 


with Bu = 0 on 3D o(x) > 0. Caugbey and O' Kelly (Reference 2-1) have shown 
that Equation (2-1) admits classical normal . modes if the following conditions 
are net: 


(1) and L,. are self-adjoint spatial operators. 

7m h 7^ 4 

(3) Tlie boundary conditions proscribed on 3D are compatible with the 
operators and L,. Under these conditions there exists a com- 

plete sec of linearly independent eigenfunctions X^(x), I c[l, «) 
such that 


/ P(x)X.X. dx = 6. . 

/ - 1 J - ij 

Jo 

I X.L.X. dx= 2w.r,.' 

/ 1 1 J - 11 

Jd 


( 2 - 2 ) 


ij 


(2-3) 



(2-4) 


f ■. 



Thus if 


u(x, t) 


y^(t)X^(x) 


(2-5) 


Equation (2-1) reduces to 


where 


^1 •" 


qi(t) 


qj(t) 



t)f(x. 


t) dx 


1 e [1, -] 



For horx)gencous initial data. 



2- Discrete Systems 
Given the system 


Mu + Du + Ku = f(t) 


u(0) = u(0) = 0 M, D, K - N N 

Caughey and 0* Kelly (Referonre 2-1) have sh.cvn that Equation 
classical norr.al Jcs iff i>.nd cc.—iuto. 


(2-7) 


( 2 - 8 ) 


(2-8) admits 



If M, D, and K are s>nffi!etrlc with H positive definite and D and K at 
least positive scaidef initc, there exists a coaplete set of ordinary eigen- 
vectors such that: 


M 6 


ij 


K 


(2-9) 

( 2 - 10 ) 

( 2 - 11 ) 


If we write 


1 C [1. Ml 


then 


u(t) = ^y(t) 


( 2 - 12 ) 


V, ■*■ v.y. •= q.(t) 

' i 1 1- i 11 1 


>'j(0) = y.(0) = 0 


(2-13) 


where 

o(t) = 4>^f(t) 



For homogeneous initial data. 
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3. Excitation of Pure Modes 

If in the case of Subsections II-A-1 and II-A-2 the forcing function is 
given by 


f(x, t) = o(x)X.(x)p(t) 


or 


f(t) = p(c) 


(2-15) 


we sec that 


qjt) = p(c)5j^j 


Hence, in the case of Subsection II-A-1, 


y.M ■ [ 

.'0 


T) P(-) dr 6 .. 


h.(t') = exp (-. r,.t) 
1 J J 


sin ^ . t 


thus 11 (x, t> 


A . ( X / V ,( 

I - i 


t;. 


(2-16) 


(2-17) 


0 


J 


i ; 
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In the case of Subsection IT-A-2, 


,j(t) - 1., 


(t - t)p(t) dT 6. 


h,(t) = exp (-U.s^t) 


sin o t 

3.. 


■j 


(2-18) 


thus u(t) = 

Hence we sec that in both cases a pure nonial node is excited. In 
particular, if p(t) = Pg cos u.t, then, as t - we have 


or 


u(x, t) = Xj^(x) 


u(t) = $ 


(i) 


cos iJt - 


J(4 - -'T ^ 


^0 ~ 
v(*’i " " f 


(2-19) 


vhore 


-1 

a. ~ tnn *“o**'Y 
^ ^ t ~ 




Using the result of Equation (2-19), wo can detorniiie and : 

4, Uonclassical Nomal Uedes 

For simplicity, we shall restrict the discussion to discrete, viscously 
damped systens. If in r.ubsection ll-A-2 m"’d and do not comnute. then 

classical norr_al moaos de not exist and it will be sho-.-n that it is ir.pes- 
slble to excite pure eigennodos by any choice of real forcing functions. It 
should he noted that even in this case it is possible to excite ^’fairly 
nun-*** i.-: ^.■nno.:..*5;. 


1 - 
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Tlu' tomuKitlon of 2N space is iis follows: 


Let 


then l't|uation (2-8) can he rewritten In the fom 


2 - Az + b(t) 


r(0) - 0 


wlie re 


aiul 



If the rat r lx A is 
tlial; 


there exists a neasluitular natiix T such 


t“^\t 


"l 



c-jn 


where 



v:-J2) 


“1 


0 
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is a diagonal nacrix of coaplex eigenvalues. The natrix T has the 
structure: 




❖ 


-S>A, 


❖ 


«A, 


where 


(2-23) 


and 




+ X.n + = 0 i c 1^1, N’j 


In general, the are connlex N vectors. The inverse of T is given by: 


where 


= 


C3 I c. 


c, . ‘ 


(2-2A) 




C, = 4. 


c, ^ c. 
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Lot 


z = Ty 


Then 


y = Ay + ^(t) 
v(0) = 0 


Q(t) = T ^b(c) 


y(t) = 


exp^A^(t - T)) 0 

0 exp^Aj^(t - t)^ 

Usinj; Kquations (2-24) and (2-25), we have 


f 


T ^b(x) di 




(2-25) 


(2-26) 


(2-27) 


x(t) = 


2Ke I exp^/.^Ct - i)) M*'^f(:) di 

Jo 


(2-28) 


/ 

/ / \ -i\r 


x(t) = 1 

Irj^. ext>/-..(t- T)j<J. j[^l 

"•(■^ ‘1* )J 


-1 


(2-29) 


Since M ^f(t) is a real vector, and exp(.‘. jt)«J> and S arc real 

r.:it rices, no ciioice of the lorcin.^ function f(t) will result in the excita- 
tion of a pi^rjL^ eigeiir:ode. It r^av also he shown that it is impossible to 
cxci'io a pair of c rplox cvuijuy.are pure eigennoios. Tlius, iinliV.e tne systen 
or Subsections Il-A-1 ami iI-A-2, nc.dal testing does not enable us to 
accurately itlenrify the eigt-nva.lues and e i genvectors of systenis wirri mm- 
c^issic:! normal r. s. I.e-pilc :ir:s fact, ii in Ib;uaticn (2-f| tfsi dam^^ing 
m:ir:x i‘ is small and the ma-r values Wt-Il separated, Kvu^lien (2-2r) can bo 


i X 1 a . . 
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Let bo such that 


0 0 


I; 


T 

& Iv^> 
0 0 







(2-30) 


If VO sot 


f(t) » cos ut 


(2-31) 


and in Kquation (2-30) 


l^ijl '-^i 


Vi, j 


(2-32) 


then after the initial transients die out. 


x(t) - 


P. cos (ut - a ) 

"LT7T""v ( 


(2-33) 


-1 "^J± 

ct = tan T 

^ u\ - (0^ 

t 


If 


yj- 


I 

. . < 
1 ] 


'} 



ird . and c avo distint-t and veil separalod, and . ‘ . . * 

‘ ‘ i . j 

tiu- nain et u et ef nnnr I i *.m 1 Jimpin.' is ta ^ 




ORIGIKAl PAGE fS 
OF POOR QUALITY 


shifts in the response vector. Each nass no longer passes through its 
equilibriuci portion at the same tine ns all the other masses, as was the 
case for classically damped systems. In particular. Equation (2-33) shows 
that if 


x(t) = 


Pq cos (tot - 




(2-34) 


Thus, for snail damping, the response is almost a pure normal node. As the 
damping increases and the separation between the eigenvalues decreases, the 
effects of nonclassical damping become stronger and the response is no longer, 
even approximately, in a pure normal mode (Reference 2-2). 

5. Useful Properties of Discrete Systems with Classical h’omal Modes 

Returning to System (2-8) in the case of Subsection lI-A-2, the Proper- 
ties (2-9), (2-10), and (2-11) can be rew^ritten in the form 


*■ [=">. ••• i®] 


= I 


r - j 


(2-35) 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Slnci* tho vectors aro .lluonrly ImJopoiulont , ttio r-»trlx <■ l;» 

nonslnjiular; theroforo: 



( 2 -) 6 ) 


U. SYSTKi KElirCTION 

Vt* shall now provo two inlorostlnj; tlioorocs. 


1. Thooron I 

A cont ir.vioiis tlyna^^ioal systors, surh as that of Subsoctlon II-A-I, Is 
givon. It oxhihits a cosploto sot I'f linearly lnilojw*n«ient viseously ilampeJ 
elassical normal nodes havinp, e 1 p.ent iinet ions X^(x), and e i p»envaliies 

I L [l, "]. Given a positive iiueper N. there exists an N^^'-orJer vis- 
Ccmsly damped discrete system exhibit ine, a complete set classical normal 
modes liav’iriit t!u* pri'p**rty th»it its i eigenvalue corresponds exactly to 

th » » 1 i 

tho i eigenvalue ot the continuous system and lurther the i eigen- 
vector corresponds to a pro]eciitni of the 1 * e igenl uiiir t i on I't the 

continuous system. That is. 




X (x > 


l. J . [1. n] X. . ‘J 


(2-57) 


‘'roi>! : 


S::;ce the c i gt'nt unc t i ons of tdu* cotK in?:ous problem are linearly inie- 





]>t’micnt , the function 
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cannot vanish identically unless a. - 0 Vi c (1, K) • Hence, there exists 
K points Xj j c (1, N) such that the vectors 


= jxj(x^)| i, j c (1. X) 


are linearly Independent. 




Since the arc linearly independent, the natrix ^ is nonsingular. Let 


•> -1 ^ -I 

M = a“(« 




where a is chosen such that 


TrM = I p(x) dx 


Using M, n, and H constructed as in Equation (2*41) to forn th.e M "order 


discrete svstcci. 


Du Ku = {(O - (t) 


u(0) = u(0) = 0 


(2*43) 
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This systca has the following properties: 

(I) There exists a cocplcte set of ordinary eigenvectors 

i c (1, N) such that if * = •*' J 


as in Equation (2-AO), then 


T 2 

^ H<fr = a I 


(2) 4> D-f = 


(3) 4>K«> = 


r ’•■■■ j"’ 

r-u-’ 


(A) 


If f(x, t) q^^(t)X^(? 

i=l 


(x)o(x) 


then q. 


Jit 


p(x)f(x, t) X.(x) d:< 


(2-4A) 


If in Equation (2-^3) u = ^y* Equation (2-->43) is reduced to- 

>'i ^‘"iS-'i ■*' Vi " 


(2-43) 


which is exactly the sane as the i^'* mode of the continuous system of Sub- 
section II-A-1 with f(x, t) given by Equation (2-AA). 


Hence, since 


;;(t) =--4>y(t) 


(2-46) 



«j(t) » 






Uj(t) £ u(x^,t) 


(2-49) 


That is, the solution of the drscretc Problen (2-38) is the projection of the 

solution of the continuous problem of Subsection Il-A-1 with f(x.t) siven bv 
Equation (2-4A). 

It should be noted that there exist, in general, infinite sets of points 
j ^ [1. N) that nay be used to define the sets of vectors 
t £ II, N). Thus there exist infinitely n.iny N''^''-order models that can be 
used to mimic the behavior of the continuous system. It is not surprising, 
therefore, that observations at N points in a continuous system do not permit 
uni.iue ident i f ical ion of the continuous svsten. 


Tlieoreia II 

Given a discrete N*^'-order dynanic.-.l system exhibiting a complete set of 
linearly independent viscously damped classical normal modes i {i. v). 

Given any positive integer N, < N. there exists an N,'“-order discrete 

dynamical system exhibiting a complete set of linearly independent classical 
normal r.odes . i r fl \-1 .i . . tb 


^1 , • liaving tlie rropercv that its 1* 


' : lut* 


™«lr t» the I ” u.o i.r*... .j.to, 


corresponds .-vact ly to a projection of 




ot liio iarj^T sysion. 


* i »;en~ 
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Proof: Since the eigenvectors of the given systen are linearly independent, 


the matrix 


.[.«>. 1®. 1®] 


( 2 - 50 ) 


is nonsingular, that is. 


1*1 i 0 


( 2 - 51 ) 


However, if this is true, there nust exist nonvanishing minors of all orders 
loss than N. In particular, there must exist at least one nonvanishing minor 


of order Let 


= S i C jl, N,j 


( 2 - 52 ) 


Wiere S is an matrix whose columns consist of cither the null vector 

0, or distinct unit vectors e., j c [l, have zero entries in all but the 
row, which has unity. The matrix 


[ ( 1 ) ( 2 ) 


( 2 - 53 ) 


is such that j^l is a minor of ^ of order * 1 ^ and so does not vanish for an 
appropriate choice of the matrix S. Let 

T 

2 - 1-1 
>L = ) 


■’2 ■ JA 

=2 = 7 


where i is irhorvn so tl)at 
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Using the M^, D^, and constructed in Equation (2-54) to fom the 
order discrete systen. 


H2U2 + D2U2 + K2H2 = f2(t) 


u_(0) = u (0) = 0 

"2 


This system has the following properties: 

(1) There exists a complete set of ordinary eigenvectors 


c [1. »,] 

such that if ^ = 



^ M2^ = a"'! 


T 



J" 

T 

*] 2 

^ ''2'^ " 



(4) If f(t) = ^ , then 




If then Equation (2-56) is reduced to: 

Zi + = q.(t) 

z^(0) = z^(0) = 0 
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Equation (2-59) is exactly the same as the i mode of the large system. 


Hence, since 


y2<t) = ^z(t) 


(2-60) 


S2i(t) 




(2-61) 


Using Equation (2-52), we have 


y2i<t) = s 


j=i 


(2-62) 


= Sx.(t) 
-1 


/. u^(t) = Sx(t) 


(2-63) 


That is, the solution of the lower-order system is the projection of the 
solution of the higher-order svsten. 


It should be noted that there exist, in general for each N, several 

nonvanishing minors. 'Tnerefore, there exist several -order models that 
minlc the behavior of the S^*^-order system. It is not surprising, therefore , 


that obser\Mtions at N’ points 


,,th , 
in a N -ord 


der discrete system do not, in 


general, permit unique identif ication of the N -order system, unless ■ N 
or there exists scr.e special structure to the system. 


Th»*orems I and IT enfi-irc* that there exist f i n i t e-d irons i on lumped param- 
eter roJois that can mimic exactly the beli^ivior of a hi glu^r-ordcr himped 
parameter or contlriuous system i;ruier anpro]>riate conditions. 


Im cr:^c tice ve eh-erv-' thv r» 


; : rs/ res:*' 


a sv'^tem at or:v a ^ j I i tiaite 


is el *' > r ! V a ^ - 1 ' t i 


mi 
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forces to the structure at only a snail finite nunber of points. By limitin- 
the number of points at which we excite the structure and observ»c the 
response, ve prejudice the outcome of any attempt to identify the structure 
uniquely. In general, if ve excite a structure and observe its response at 
points, ve can uniquely identify only an -order lumped parameter 

model. 


EFFECT OF MODELING ERRORS ON PREDICTED RESPONSE 
Given the System (2-6^), 


f f(t) t ^ T 

b: + Dx + Kx « / 

( 0 t > T 


x(0) = i(0) = 0 


It is sometimes more convenient to write the equations in Ist-order 


= Az + 


|s(t) t 


0 t > T 


z(0) = 0 


■fH: -:i 


i " ) 
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Suppose now that ve have an analytical nodoi of System (2-65). 


j; ■ "s * 


I h(c) t ^ T 

( 0 t > T 


( 2 - 66 ) 


y(0) = 0 


where 


B ^ A h(t) =• s(t) 


(2-67) 


U'e wish to know whnt errors are induced in the solution by Kodeling 
errors in A and g(t). Let 


(2-68) 


thus 


thus 


dw 

^ = Aw + (B - x\)v + js(t) - h(t)[ 
w(0) = 0 


(2-69) 

(2-70) 


•Kin(t,T) 

•(c) = I exp A(t - t) j^(3 - A)v(t) + ^r(t) - h(i)Jj d: (2- 


70) 


tnus 


•-yinft,-: ) 


|w(t),i 2 


c?:p A(t - :) 




A) : 


- h(') ; 


’ 1 ) 


0 
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^"ow 


Ilexp (At) I I - exp(-aj^t) 

Ijexp (Bt)ll ^ expC-a^t) a, > 0 

Let 


Then 


i|w(t)!| ^ 



M = Max (Mj, M 2 ) 


i = Min (a^, 


d = Sup I |g(t) - h(t) 1 1 


t 

Max 


X ISupi jg(c) I 1 , Sup llh(t)il 

L t t 


a! t Sup I |y(t) I 1 + 


dj '] j^exp (-a 


(t - t ) 


vhere 


Nov 


* 

t = t if t < T 
= T if t - T 



Min(t,T) 

exp 





d' 
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^ Sup I |h(T)i I |t-xp ^-n(t 

liw(t)l| - L ^cxp ^-a(c - t*)^ 

+ ^ d ^oxp |-a(t - t )j - fsp 

There are two cases of special interest, 

1, Case (a) - Short Transient toads 
If aT << X, then Equation (2-7S) yields: 

siipi|w(t)i| - (>rr)‘kllB - a|1 + mtj (2-79) 

t 

If MT 0(1), then the two terms d and k ] |b - a| j are of equ*il inpertanco. 
!U-nce, the errors ]|b - Aj 1 and |[^(t) - h(t)jj are of equal sie^n i i icance 
and tlie system parameters need not be known with anv i:ii;her aceuracy than 
t he f or e i up, f un c t i ons . 

2. Case (b) - Tersistcnt Excitation 

If ^T 1, then F'quation (2-78) yields: 

Supi!w(Ol| k Mb - a|: d (2-80) 

Since in many stnictur.il dynamics problems the dampine, is small, the 
LiTiTi M/a boccr.u‘s v*‘rv Uirco eoT?.p,s r^d to unity. In thi:» case, t!ic first tern, 
k \\\ - A. I assumes much c;ieater i rticrt ance than the St-coiul tt rm d; this si.ovs 
r!ie pv>ssihle effects of resonance- Tluis ve see that u:u!t*r ]'‘.rsis:ent o.cit.i- 
r:o!i, svst.-m errv'rs <'an r!av a domln.int roll; hen.e, the sy-^:i ;>.\rer^:o:s 
m-:st h,- J.-: ;:u'C villi a ruv'h hii^her .:ccuracv !1 mu : ho nucin.: ' uiu' ; : c!Vs . It 
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Is interesting to note that Chen and Uada (Reference 2-3) established a 
sinilar result using pcrturbational analyses. 


D. EFFECTS OF NEGLECTING HIGHER MODES 

Let us tirst consider the case of an order discrete System (2-8) in 
which the forces f(t) are basically low frequency in nature 


Mu + Du + Ku = f(t) 


u(0) = u(0) = 0 


where and M connutc. Then 




u(t) ^i<t) 


r ^ exp (-l: , r, . ( r - r) ) 


sin (t - t) q^(T) di 


where 


(2-8A) 


q(t) = ^ f (t) 






(2-85) 


"i *i ‘ ' ’i 


c (1, :.•) 


■rr-c?' 
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If the frequency spcctrun of f(t) contains only frequencies which arc 
snail conparcd to all but the lowest eigenvalue then 


for i > N*, N* < N 


In this case it is reasonable to approxinate the solution u(t) by the 
truncated series. 


( 2 - 86 ) 


u(t) = u^(t) 


(t) 


i=l 


(2-S7) 


If ve look at a typical term in the remainder of the series. 


Jo 


,(t - t) q^(T) dt 


■/ "j<'> ’j 


(t - l) dl 


(2-S8) 


Expanding - t) in a Taylor's scries about t = 0, we find 

2 (t) = q (t) / h (t) d- - q (C) / t!i.(t) dT 

Jo ' ' Jo ' 


q,(t) 2 

+ / T hj (t) dT + - • - 


(2-89) 
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If > 1, then 


hj(T) dT = 


/ Thj(T) dT = ~f 

h ‘"j 


f 2 1 - 

/ Th.(t)dT-— ^ 

lo “i 


q.(t) 2,^ 

< gax 




i.o., q^(t) is lov freq’jency cor:pared to , then 


j c + i, 


(2-93) 
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In this case ve iiay inprove the approxicate solution of Equation (2-87) 
by adding the additional terns of Equation (2-93). 


u(t) ^ 


(2-9A) 


(s) 

4 (t) = -^ i c (1, N) 


(2-95) 


Equation (2-94) nay be rewritten. 


j=l L J J 


(2-96) 


(t) _ *T^. . 
q = * f(t) 


(2-97) 




(2-9B) 


2 


y.l:- 



(2-99) 

k~ 

(2-100) 


(2-101) 





where u ‘ (t) is the 

f(t). 


"static** response of the system to the applied forces 


If we write 

* 

N 

(2-102) 

1=1 

then 

u(t) = +u^®^(t) (2-103) 

Thus we see that the total response is approximated by two separate 
(d) 

terns; u (t) consists of the d>Tiarnic response of the active modes less the 
’’static" response in these inodes and u^^\t) the "static** response of the 
whole structure to the applied forces f(t). These results are identical to 
what ir sometimes called the node acceleration method. 


As a measure of the error in neglecting the higher modes we have 


I ly(t) 1 j - I |u_(t) i 1 

!iu(t)| I 


(2-104) 


This error tends to be smaller for the case of persistent excitation where 
resonance nay occur, and higher for short transients where resonance has no 
chance of occurring. 


E. ERRORS IN EIGENVALUES /\ND EIGENITCTORS 

From an analytical point of view, errors in eigenvalues and eigenvectors 
arc usually the result of codding errors or the use of too crude a level of 
(iiscrer iration of the continuous system. In practice, it is usually possible 
to model the lover eigenvalues and eigenvectors of aerospace structures with 
an accuracy of five to ten percent, or better. The eigenvalues are usually 
noro accurately modeled than the ci'/f:r.vectorr^ though this may 
a prcblvm rolatt'd to the difficuirios encountered in testing and measuring 
o 1 cer.va _ues and e i gvnvt.o: tors. liiis 'Uhiec.t :;c treated in ro'-rc* detail 

in Section III of this report. 
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F. ERRORS IN MODAL "FORCE" COEFFICIE.TS 

Errors in eigenvalues and eigenvectors have a significant effect on the 
modal "force” coefficients; however, even if the eigenvalues and eigenvectors 
of the first N nodes are Vaiovn exactly, errors in the modal force coefficients 
will still arise due to discretization effects. For example, a central dif- 
ference approximation may be used to calculate the curvature of a beam using 
the discrete displacements of the bean. 

2 , 2 

Thus 3 u/3x is approxi:iated by 


u + u , - 2u 
n+1 n-1 n 


(2-105) 


where u = u(x), x =nh, and h is tlie mesh s^sacing. Kow 
n n n . 


Thus 


If 


thcr. 


ft ■^11 — ^ 7 

n+1 n-I ^ n 3‘u . h 3^u . ^,.4, 
^ = — ^ + — - — — + 0(h ) 

h (?x dx 


(2-106) 


1 

n+1 n-] 

2 

^_2„n-h —2 

.X 

^ 2 

< " - 

/ 

^ *♦ 

1 

,2 

rA 


“ 12 

.2 
0 u 





dx^ 


max 


max 


max 


(2-107) 


u(x) = A sin 


in — ^ 


(2-103) 


I n-1 


2u - h ;-l 

n 2 ^ 

%:< I 


2 


I % 11 1 


(2-10'^ 


\:y. I 
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Now X/h =* the number of mesh points per wavelength; 


^ -2 

u + u . - 2u - h 
n+1 n-1 n ^2 

<x 


Is.' 


nax 


3N 


(2-110) 


m 


max 


since 



( 2 - 111 ) 


The relativ 2 error in the bending stress is also given by Equation (2-110), 

Thus, the relative error in the bending stress increases rapidly as , the 

m 

number of mesh points/wavelength, is decreased. Since we have shown in 
Section II-B that it is theoretically possible to construct discrete models 
whose eigenvalues agree exactly with the first li eigenvalues of the continu- 
ous system, and that eigenvectors are projections of the eigenfunctions of 
the continuous system, it should not be too surprising that the accuracy of 
modal “force” coefficients is usually lever than that of either the eigen- 
values or the eigenvectors. 

In addition to the discretization errors discussed above, additional 
errors arise because of nodal spill-over caused by experimental difficulties 
in obtaining pure icodal excitation, and the fact that the real structure nay 
not admit classical normal modes; these errors will be discussed in more 
detail in Section III of this report. 





SECTION III 
DYNAMIC TESTING 

Though it is possible, using nodem analytical and numerical techniques, 
to model the dynamic behavior of a structure to any desired degree of 
accuracy, few engineers would be happy to commission a new space vehicle 
without at least some limited dsTiaaic testing. The main reasons for this 
conserv'ativc approach are: (1) while it is possible in theory to model the 

structure accurately, it is usually too costly or too tine consuming to do 
so; (2) it is ver>- easy to omit sone significant effects such as geometric 
or material nonllnearities in modeling the system. For these reasons, most 
engineers require at least a limited program of dx-naaic testing to **qualify** 
the analytical model- If the limited program of testing yields good agree- 
ment between measured and predicted values, the engineer is happy; unfor- 
tunately, if the agreement is poor, as it frequently is, the engineer is left 
In a quandry. It has been suggested that the test data be used to update 
the analytical model and so increase the precision of the analytical predic- 
tions. This is a very useful technique and can yield good results if prop- 
erly applied. First of all, it must be pointed out that due to the non- 
uniqueness of the identification process, it cannot be used to identify the 
parameters of the complete structural system. It can, however, be used to 
obtain updated estimates of the natural frequencies, damping factors, and 
mode shapes of the finite numbers of nodes observed, and hence to make better 
estimates of the contributions of these modes to the response of the struc- 
ture. While this approach is useful in improving the analytical and predic- 
tive capabilities for a given physical structure, it cannot help improve the 
analytical and predictive capabilities for new and unbuilt structures. 


MODAL TESTING 

As shown in Subsection II-A-3, discrete and continuouj; systems exhibit- 
cl*'issicai norr.ji -odes are capable of being excited in pure normal modes 
shu-Ti in Subs(-ction TI-A-4, discrete systems exhibiting nonclassicai nor— 
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and the eigenvalues well separated, the lower nodes of any viscously danped 
Structure can be excited in rather good approxitaations to pure nodes* 

In structural dynanics, it is cocrsaonly assumed that the system possesses 
classical normal nodes. While this is seldom strictly true, the damping in 
many aerospace structures is often quite small, and so fairly pure modes can 
be excited, particularly in the lower modes. 


Consider the N -order discrete system that exhibits classical normal 


modes: 


Mu + Du + Ku » f(t) 


u(tQ) * “ Q 


f(t) = C cos wt, tQ > 


Since Equation (3-1) has classical normal modes, the solution of Equation 
(3-1) with Equation (3-2) can be written 


u(t) =• exp (jc-it)j 


‘^l^cre is the nodal natrix (2-12) and R Is the response nstrix: 


exp (-ja.) 


J(4 - 4 


4- 


’vsjcrc * Is defined in F.cuation (2-19). 


(3-5) 




Then 
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u(t) = Re ^<frRq exp 


(3-6) 


If 


3 = %£ 


k-k 


(3-7) 


then 


u(t) = cos (- :t - ak) 



'4 - ^ 


(3-8) 


Condition (3-7) requires that 


C - 


(3-9) 


However, d> is unknown at the beginning of the test. 

We observe that if - 5 r.^ « Vi and the eigenvalues are distinct 

and well separated, then \k^\ » i * k; thus; 


u(t) = 


cos (-t - ^J_) 


(3-10) 




This suggests that an iterative scheme can be used to generate pure nodes 
Using Equation (2-6), we can define an iterative schene. 


uj^(t) = ?e Uj^(t) 

u, (t) = 4-?. q. , «‘>;p vj.;t) 
c (t) = " V constant 


(3-il) 


' 12 ) 
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Usir.^ Che properties of classical noraal modes, we have 


I 


9k ” ^3k-l 


where 


It exp (-.i'lj) 


2 2 \'‘ 2 


Kj. =0, I t i. 


3i = ^3o 


32 = = ^“9o 


3„ = = K 2o 




1 Vi 5^ 


;>! i: V.;- select • ex') (-j ), ) s’.ic!: Chat 


thon 


PC>C^ Q'-’^'*-’ • ' 


lia R-‘ = E 

jj.«= 


(3-19) 


IChcre E. are nonnejative, definite natrices satisfy the following property: 


"•i^j = "i'ij 


(3-20) 


Thus, using Equation (3-16), wc have 


° " '’oi^i 


(3-21) 


Using Equations (3-21) and (3-11), we find 


lia u (t) = lin*Rq,_ exp (j^^t) 
k 


cy.ry [ j (u:t a ) J 

_ ,(i) L_ 




(3-22) 


lin 11. (i) = lin ii(t) 


-k 


k • 


(i) cos (-:t - 


2 2 
"i 


)t 


(3-23) 


Thus* -hi iterative process converges to .1 solution ])roport :• nal to the pure 
--he iron vhich the natural frecuenoy , ani the p-irancter 

i) 


- i V 


e-h:,iini^h in additii^n to ‘ 



It should be noted in passing that what has been said for the discrete 
System (2-8) is also true of the continuous System (2-1). Any system, con- 
tinuous or disoretc, that exhibits classical nonnal codes lends itself to an 
iterative scheme, like tliat above, which converges to a pure node. 

B. ERROI^S IX MODAL TESTING 

1. Nonclassical Noma I Modes 

Even if the system under test can be regarded as an N-degree-of-f reedom 
discrete system it is highly unlikely, in practice, that the system will 
possess classical normal modes. Tlius, theoretically, there exists no choice 
ot toicing functions that excite pure nonnal codes. As previously pointed 
out. If the system damping is snill and the eigenvalues well separated, 
relatively pure modes of vibration may be excited. If the system damping is 
iiv»t small and the eigenvalues closely spaced, as often happens in real struc- 
nires, it nay be impossible in practice to excite even relatively pure modes 
of vibrations in all but tiu? very lowest nodes. 

2. Impure l!cdal Excitation 

Even if tlie system under tost can he regard«*d as an N-degreo-of-f reedom 
discrete systt'Cj possessing classical nom:il moties, excitation of a pure code 
ot vibration rec;uires that eacl: mass In the system be excited hv a force 
proportional to that mass and to the nodal displacement of that mass. In 
Equation (3-1^, we have shouni that an iterative tecJinique can he used to 
achieve this end, provided the mass matrix of the svstem is known and pro- 
vided that we have the means to apply forces to each mass. Wiiile we will 
seldom know tne mass matrix precisely, we often have .uK‘quately good esti- 
mates; however, wo seldom have N-fnrcc transducers available to conduct the 
test. In the case of a continuous system, modal excitation sljould also be 
continuous, and at least an adocpiate <Hscrete approx imat ion to a com inuous 
d i St r i hut ion o: forces an even more d i t 1 icu 1 1 t ask to acci^mp I i s!i . 


Mca-urv m^.nt Errors 
Assuming, as in Snhscctions 
rt‘e o: I rci is a gOv-d n:v>del of the 

‘‘Itn c! T.car'urinL; d i >.;> 1 acemcTU r, a 

;i c: tho } ore in;- i Muction, Vi 


ai]d ni-U“2, tliit the discrete 
s V s t L-m , 1 1 ic r e still r ema i n s t ] k* 
'• ; Jir:ts c-.e i . 

St C'! nil, it i liot .isa.t 1 : n 


,-li 



N displaocacnt transducors available for a lost; true. It Is possible to^uso 
S ^ K tratisduoors. and novo then around the structure, but this is a t i=o 
cLsuniu, procedure and preatlv increases the cost of the test. However, even 
if vv have N transducers available, and even if ve could excite pure normal 
modes, there still exists the c,ucst ion of n<>as«ro=ent error, particularly 
errors in phase ceasurenonts, which are notoriously difficult to make with 

accuracy* 

4. Fffocts of Discretization or Condensat ion 

Aerospace structures are alnost always continuous in nature, or, at best, 
very-hij;h-order discrete systems; however, for purposes of analysis and test- 
inc, we rust discretize the conliuuous structure or condense the high-order 
discrete svstem to obtain a n,inageable system, in Subsection III-B it was 
shown tKit if the original system, continuous or discrete, exhibited classical 
normal m,xles, it was possible to construct an N-degree discrete model whoso 
oicenvalne coincided exactly with the first N eigenvalues of the original 
svstem and whose eigenvectors were projections of the eigenfunctions or 
oigenvoctors of the original system. It was further shown that tor certain 
cl.ssos of excitation, the response of the model exactly mimicked that of the 
original structure. Despite these very useful properties, it should be clear 
that the model is not one to one with fne original structure. This tact 
shews up immediately in modal testing. hot us suppose that we have a con- 
tinuous structure and that we assign to it N coordinates X. i = 
that we shall nuike measurements and apply forces only at these X points. 

In E.iuation (2-^3). let f(l) be given bv 


I j u) 


fit) = 


( 3 - 24 ) 
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The points of application of these forces correspond to the points 
i c (1, N') of the continuous systca of Equation (2-1). Thus 

N 

f(x,t) = ^ ?. (t) 5(y. - X.) 

i=l 


q,(t) = 


I 5j^(x)f(x,t) dx 

.h ~ ~ 


= Vx.(x.) f. 

-J J 


(t) 


j=l 


With Equations (2-39) and (2-40), if i c (1, N) , then 



bat 


f'(t) = >'$q(t) 


Thus 

o(t) = 4>^f(t) = $^M$q(t) 


q(t) = g(t) 


Drovi.j»fJ ?h::l 


(3-25) 


(3-26) 


(3-27) 


(3-28) 


(3-29) 


( 3 - 30 ) 


(1 



Ke note that 
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a,(t) 

j=l 


(3-31) 


1 c IN + 1, “J 


In general. 


q.(t) i 0 


(3-32) 

(3-32) 


For exanplc, if 


cj(t) = qj^(t)e|^ 


k e [l, II ] 


(3-33) 


Thus q.(t) = ^ 


q,(t) = > :;,(x.)f.(t) 

^ Z— ^ ^ -j j 




^ 0 in general 


1 > N 


(3-34) 


th 

Thus, for the- first U nodes, onl^*” the 1: node is excited. However, 

there exist higher-order nodes, t > N\ which are excited. If k is much 

smaller than N, then the separation in eigenvalues is usually such that the 

response of these ’haliased * nodes is sniall cfcip.ired to the response of the 
th 

rr.ece, particularly Lf the 1 recueucy of e>:cItation is close to the natural 

- _ t h 

irequency of the k rode. If, lU’wever, the iroquency of excitation is close 
to that of one of the "aliased” nceies, serious errors can result. 



c, OTHER ident:ificatio:; techkiques 

Since nodal testing as an identification technique is restricted to 
classically danped systens, or at least systCRS with small damping, it is 
reasonable to ask if there are ocher identification techniques that could 
do a better job. There exists a variety of identification techniques, both 
parameter and nonparameter ; however, if one wishes to identify the mass, 
stiffness, and damping natrices, one is faced with a fundamental linitation: 
Che number of points at which neasurements are made must, in general, be 
equal to the number of degrees of freedom of the structure. Unless this 

is done, the solutions obtained are not unique. Since, as already pointed 
out, aerospace structures are almost always continuous, or at least have a 
large number of degrees of freedom, unique identification of c!ie structure 
is virtually impossible. One has to conclude therefore that, at least for 
structures with small damping, modal testing is probably as good an identi-- 
fication technique as any available. 



SKCTION !V 

CONXLUSIOSS ANl) RECO>CIl'.\UATIONS 


Tlie object of this report h.is been to cxanlne the problcas of analyses 
and t*-sting of aerospace structures and the difficulties of correlating 
the results. 

A. CONCLUSIONS 

The following conclusions have emerged fnra tills study; 

(1) Mv>dal testing (and indeed any other teclinique) cannot be used 
•IS a aethod for unifjiiely deternining t!ie mass, damping, •ind 
stiffness n:itr{ces of real aerospace structures, which are 
usually continuous in nature, and the nodal testing thereby cannot 
provide a means of improving the analytical techniques for deter- 
mining dynamics response. Modal testing is an extren^'lv useful 
tool for obtaining accurate measures of the eigenvalues and eigen- 
vectors of tlie lower modes of tl;e structure. These measures can be 
correlated with analytical results, or provide the basis for dis- 
crete models of tlie structure, which nay be extremely useful in the 
analysis of stability a:>.! control. 

(2) Analytical techniques are capable of modeling dynamic structures to 

any desired degree of accuracy. It is clear that eigenvalues can 
be pre\i feted with a higher degree of accuracy than can eigen- 
vectors, and that eigenvectors can be predicted with a higher 
degree of accuracy tfian that accuracy with which the rx^dal force 
CiH'fficient can be predicted- Accurate analytical prediction of 
^'t^esse^• an.! lorces will rt'ju:r<* a finer level 4 >f d lsci‘i‘ t i cat it>n 
than will the predictiiui of eigenvalues and t^igenvectors. li accu- 
r.iii' araivtical ; i m t ien*- ».m ? tres> ar;d i\>rv es are recviin’d, tlie 

additioiuil 4'tfer: .uul exrer.se of using liner and finer r.e^iH-s must 
be ar.'epled. 


(3) Accuracy of nodolin>; is a central question in all mathematical 
model iuj;;: "tUven that tlie data luis only finite accuracy, lu^ 

acciirate need t!ie node! he to obtain acceptable accuraev i:i tlie 
response?” Tins topic was dealt with at some length In Sub- 
section 11 -C, and the ansver depends on the fora v>f excitation. 

For short transients, the nuxlel need not be any more accurate than 
the input data; for persistent inputs that create the possibility 
of resonance, the model must he specified with a much higher 
degree of accuracy than the input data. These results are consis- 
tent with the results obtained by Tht'n and Wada (Keforence 2-3). 

KFC0MM12;iUTI0!%*> 

As a result of the present study, two recommendat ions emerge: 

(n To determine stresses and forces In aerospace systems, an appro- 
priate level of d iscret i^at ion mist ho used, even If this is much 
finer than would he uscvi for determining eigenvalues and 
eigenvectors. 


(2) While it is virtually impossible to "idei^tifv” the structure 
uniquely from the result of moaal testing, such tests provide 
a valuable clieck on the analytical method and can be used to 
provlvle an accurate, discrete model of thv'-svstom for use in 
studit's 02 Stability and contritl. Modal testing lias an addi- 
tional virtue that should be exploited to the fullest. There are 
many physical phenomena, such as the sloshing of fuel in a spinning 
spacecraft that are rather difficult to model accurateiv since 
boundary layer friction and dissipation play a central role. In 
such a l ase, modal tests nt t lu* phvsical system can east I v provide 
the data on which to base an analytical model. This modal testing 
tlion beco.v s part and parcel oi the modeling techni(jt:e in wdiich 
some parts oi t l:e stricture are nodeKJ ab initio, and some parts 
tn\.>de!t“d oTi the basis oi t hi* m 


I'd* 1 t s t . 
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AN EXA.MPI.E 




Consider the problem 
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9^u 

9t^ 


+ 26 ~ » —• + f(x, t) 0 < X < 1 


3x 


u(0, t) = u(l, t) = 0 


(A-1) 


u(x, 0) = u(x, 0) » 0 


The eigenvalues and eigenfunctions of Equation (A-1) 


arc 


= (ir)^ i c [ i. "') 


X^^^(x) = /f sin i-sx 


(A-2) 


If ve write 




1=1 


(A-3) 


then 


Zj + 26Z^ + ^“z, = q. (t) 


X 1 


(A-A) 


where 


Uj(0) = ».(0) = 0 


and 
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Discretization of Equation (A-1) 

If we use central difference spatial discretization or constant mass 
matrix finite elenent spatial dlscretizat ion. Equation (A-1) beco^ics 

Uj + 26uj^ + (X)^ “ l'j(t) (A-6) 


where 


u^(t) " u(lh, t) ?^(t) “ f(lh, t) 
h = • 1 

The eigenvalues and eigenvectors of Equation (A-6) are 


(A-7) 


A ,A , 2 i- 

- 4:. sin 


.(1) . /2 ijz 




(A-8) 


Let 


Ji) 


,(i) 

1 . 

>(i) 


AO 

•^X-l 






(A-9) 




I 



Let 
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u = ❖ Z 


where 



(A-10) 


2 


+ 2£Z'j + n“ Z^ = qj,(t) 


(A-11) 



= $t'(t) 


j = l 


Conpnrisons o: Equations (A-3) and (A-9) shovp that 


=<» . x“> (i; 


(A-12) 


2 :: Sin ^ 
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r, 

i. ■ 


-] 


sin ^ 


(A-13) 







OF poen Q'JALITY 

Thus, in this case, the eigenvectors of the approxinating system are 
projections of the first N eigenfunctions of the continuous system. 


For the sake of illustration, let us use H ■» A. Equation (A-13) then 

gives 




(A-14) 


From Equation (A-14) 


= 0.0255 

* 0.0997 

* 0.2158 


1 (A-15) 


This shows clearly how the errors Irtrease with node order. If Equation (A-h) 
is vritten in natrix forn. 



(A-16) 


(A-17) 



ORIGINAL PAGS !3 

OF POOR Q'JALrrV 






Let us now use Theorem 1 to construct a 
same eigenvalues of the continuous system. 


I’^^-cn.er system having the 
Let US select 


Xj “i i «:[l. 3] 


Thus 



(A-18) 


(A-19) 


(A-20) 


(A-2I) 



OF POOR QUALITY 


In this case we have chosen a“ so that M is the sane as that obtained by finite 
difference. Thus, 


and 


11.103 


-6,9785 
- 1.2337 


-6.9785 

12.3370 

-6.9735 


1.2337 
-6.9785 
11.103 . 



(A-22) 


(A-23) 


Coeparison of Equations (A-17) and (A-21) shov that both rocels have 
the saise M and D natrices (they were so constructed); however, the iC natrices 
of Equations (A-23) and (A-27) shov ccnsiderable differences: 


(1) is a full matrix, while K is a tridiagonal Jacobi matrix. 

(2) The nagnituaes of the elements of the two K matrices markedly 
different. 


Thus, ve see that if modal testing is used to identify a 3^^-order model 
of the continuous system, the natrices of Equation (A-17) would result. If the 
stiffness matrix is compared with K, the matrix obtained from finite dif- 
ference or finite element d iscret ization, ve see that they're r.ot even close. 
Hence, ve see that o stiffness matrix obtained from nodal testing cannet be 
used to check that obtained from systematic analytic reduction techniques 
such as finite difference or finite element. 


OF poor. QUALITY 

To illustrate aliasing, suppose that g(t) in Equation (A-lh) is given 


by 



(A-2A) 


vhich will excite only the first ncnlc of oscillation of Equation (A-16). 


Using Equations (3-25), (3-26), and (3-27), we have 




^.(t) = •'2 I sin ^ /2 sin + si 


2i*r 


31 


cos v't 


sin 


T [• “ V 


COS U't 


(A-25) 


(A-26) 


Thus 


q^,(t) ^ 0 unless i = 8k - 1 


(A-27) 


k = 0, 1, 2 ••• 


qj,j.,l(t) = (±1) -1 cos *t 


(A-2S) 


Thus 


qj(t), q^ft). qg(t) etc. 0 


If 0 - J 4*, only the first mode will be stroi’.gly excited. If, hevever, 
0 - 8". the and 7^'* node of the continuous structure can be 

strongly excited. 


For structures with veil separated eigenvalues and excitation 
restricted to the bandwidth of the first N inodes, aliasing does not present 
a serious problen. Some structures, such as shell-like structures, tend to 
have rather closely spaced eigenvalues and, in this case, aliasing becomes 
a more serious problem. While this discussion was restricted for simplicity 
to the case N = 4, the sane features show up for all values of N. 








